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Hadronic r decays provide a clean laboratory for the precise study of quantum chromodynamics (QCD). 
Observables based on the spectral functions of hadronic r decays can be related to QCD quark-level calculations 
to determine fundamental quantities like the strong coupling constant, quark and gluon condensates. Using the 
ALEPH spectral functions and branching ratios, complemented by some other available measurements, and a 
revisited analysis of the theoretical framework, the value a s (mj) = 0.345 ± 0.004 cxp ± 0.009th is obtained. Taken 
together with the determination of cts(M% ) from the global electroweak fit, this result leads to the most accurate 
test of asymptotic freedom: the value of the logarithmic slope of qJ 1 (s) is found to agree with QCD at a precision 
of 4%. The value of a s {M%) obtained from r decays is a s (M%) = 0.1215 ± 0.0004 cxp ± 0.0010 th ± 0.0005 cvo i = 
0.1215 ± 0.0012 . 



1. Introduction 

The r is the only lepton of the three-generation 
Standard Model (SM) that is heavy enough to de- 
cay into hadrons. It is therefore an ideal labora- 
tory for studying the charged weak hadronic cur- 
rents and QCD. Observables based on the spec- 
tral functions of hadronic r decays can be re- 
lated to QCD quark-level calculations to deter- 
mine fundamental quantities like the strong cou- 
pling constant, quark and gluon condensates. We 
report here the results of a QCD analysis of the 
final ALEPH spectral functions [1] using a revis- 
ited theoretical framework [2]. 



2. Tau hadronic spectral functions 

2.1. Definitions 

The spectral function v\ (a\, do), where the 
subscript refers to the spin J of the hadronic sys- 
tem, is defined for a nonstrange (|AS| = 0) or 
strange (|AS| = 1) vector (axial- vector) hadronic 
decay t~ — > V~v T (A~v T ). The spectral function 
is obtained from the normalized invariant mass- 
squared distribution (l/N v / A )(dN v / A /ds) for a 
given hadronic mass ^/s multiplied by the appro- 



priate kinematic factor 

Til?' 

vi(s)/ai(s) = . . , ? T t 



B(t~ 



V-/A- v T ) 



V/A 



N V /A ds 



1 - 



a (s) = 
dN A 



EW B(t 


— > e 


m T J \ 


2s 




B{t- - 





v e y T ) 
-i -l 



(1) 



6 | | 2 Sew B(t~ 



V e V T ) 



1 



N A ds V~ m 2 T J ' (2) 
where Sew accounts for electroweak radiative 
corrections [3]. Since CVC is a very good ap- 
proximation for the nonstrange sector, the J = 
contribution to the nonstrange vector spectral 
function is put to zero, while the main contribu- 
tions to ao are from the pion or kaon poles, with 
(l/N A )dN A /ds = S(s - ml K ). They are con- 
nected through partial conservation of the axial- 
vector current (PCAC) to the corresponding de- 
cay constants, fn,K- The spectral functions are 
normalized by the ratio of the vector/axial- vector 
branching fraction B(t~ — > V~ /A~v T ) to the 
branching fraction of the massless leptonic, i.e., 
electron, channel. The direct value for B e and 
the two derived values from £> M and r T using lep- 
ton universality are in good agreement with each 
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other, providing a consistent and precise com- 
bined ALEPH result for the electronic branching 
fraction, 

= (17.818 ±0.032)% . (3) 

Using unitarity and analyticity, the spectral 
functions are connected to the imaginary part of 
the two-point hadronic vacuum polarization func- 
tions 

n^ = (-^V + ^)n^(, 2 ) 

+q^U%(q 2 ) (4) 

of vector (U^ = V? = q^qi) or axial- vector 
(Uij = A?j = ^7^75^) color-singlet quark cur- 
rents, and for time-like momenta-squared q 2 > 0. 
Lorentz decomposition is used to separate the cor- 
relation function into its J = 1 and J = parts. 
The polarization functions II^ 1/ [7 (s) have a branch 
cut along the real axis in the complex s = q 2 
plane. Their imaginary parts give the spectral 
functions defined in (1), for nonstrange quark cur- 
rents 

lmU ud,v/A( s ) = ^i/«i( s ) , 

Imng A ( S ) = ^ao(s). (5) 

The analytic vacuum polarization function 
^-ijjj(q 2 ) obeys, up to subtractions, the disper- 
sion relation 



where the unknown but in general irrelevant sub- 
traction constants can be removed by taking the 
derivative of Uij t u(q 2 ). The dispersion relation 
allows one to connect the experimentally accessi- 
ble spectral functions to the correlation functions 
Il^ u (q 2 ), which can be derived from QCD. 

2.2. Inclusive nonstrange spectral func- 
tions 

2.2.1. Vector and axial-vector spectral 
functions 

The inclusive r vector and axial- vector spectral 
functions are shown in the upper and lower plots 



of Fig. 1, respectively. The left hand plots give 
the ALEPH results [5-7] together with its most 
important exclusive contributions, and the right 
hand plots compare ALEPH with OPAL [4]. The 
agreement between the experiments is satisfying. 

The curves in the left hand plots of Fig. 1 rep- 
resent the parton model prediction (dotted) and 
the massless perturbative QCD prediction (solid), 
assuming the relevant physics to be governed by 
short distances. The difference between the two 
curves is due to higher order terms in the strong 
coupling (a s (s)/Tr) n with n = 1, 2, 3. At high en- 
ergies the spectral functions are assumed to be 
dominated by continuum production, which lo- 
cally agrees with perturbative QCD. This asymp- 
totic region is not yet reached at s = m 2 for the 
vector and axial-vector spectral functions. 

2.2.2. Inclusive V ± A spectral functions 

For the total v\ +a\ hadronic spectral function 
it is not necessary to experimentally distinguish 
whether a given event belongs to one or the other 
current. The one, two and three-pion final states 
dominate and their exclusive measurements are 
added with proper accounting for anticorrelations 
due to the feedthrough. The remaining contribut- 
ing topologies are treated inclusively, i.e., with- 
out separation of the vector and axial-vector de- 
cay modes. This reduces the statistical uncer- 
tainty. The effect of the feed-through between r 
final states on the invariant mass spectrum is de- 
scribed by the Monte Carlo simulation and res- 
olution effects are corrected by data unfolding. 
In this procedure the simulated mass distribu- 
tions are iteratively corrected using the exclusive 
vector/axial- vector unfolded mass spectra. Also, 
one does not have to separate the vector/axial- 
vector currents of the KK-jt and KKtttt modes. 
The v\ + cti spectral functions for ALEPH and 
OPAL are plotted in the left hand plot of Fig. 2. 
The improvement in precision when comparing to 
a sum of the two parts in Fig. 1 is significant at 
higher mass-squared values. 

One nicely identifies the oscillating behavior of 
the spectral function and it is interesting to ob- 
serve that, unlike the vector/axial- vector spec- 
tral functions, it does approximately reach the 
asymptotic limit predicted by perturbative QCD 
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Figure 1. Left hand plots: the inclusive vector (upper) and axial- vector (lower) spectral functions as 
measured in [5]. The shaded areas indicate the contributing exclusive r decay channels. The curves show 
the predictions from the parton model (dotted) and from massless perturbative QCD using a s (M§) = 
0.120 (solid). Right hand plots: comparison of the inclusive vector (upper) and axial- vector (lower) 
spectral functions obtained by ALEPH and OPAL [4]. 
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Figure 2. Inclusive vector plus axial-vector (left) and vector minus axial-vector spectral function (right) as 
measured in [5] (dots with errors bars) and [4] (shaded one standard deviation errors). The lines show the 
predictions from the parton model (dotted) and from massless perturbative QCD using a s (-M§) = 0.120 
(solid). They cancel to all orders in the difference. 



at s — > m^. Also, the V + A spectral func- 
tion, including the pion pole, exhibits the features 
expected from global quark-hadron duality: de- 
spite the huge oscillations due to the prominent 
7r, p(770), ai and p(1450) resonances, the spectral 
function qualitatively averages out to the quark 
contribution from perturbative QCD. 

In the case of the v\ — a\ spectral function, un- 
certainties on the V/A separation are reinforced 
due to their relative anticorrelation. Similarly, 
the presence of anticorrelations in the branch- 
ing fractions between r final states with adja- 
cent numbers of pions increase the errors. The 
Vi — ai spectral functions for ALEPH and OPAL 
are shown in the right hand plot of Fig. 2. The 
oscillating behavior of the respective v\ and a\ 
spectral functions is emphasized and the asymp- 
totic regime is not reached at s = m^. How- 
ever again, the strong oscillation generated by the 
hadron resonances to a large part averages out to 
zero, as predicted by perturbative QCD. 



3. HADRONIC TAU DECAYS AND 
QCD 

3.1. Generalities 

Proposed tests of QCD at the r mass scale [8- 
11] and the precise measurement of the strong 
coupling constant a s , carried out for the first 
time by the ALEPH [12] and CLEO [13] collabo- 
rations have triggered many theoretical develop- 
ments. They concern primarily the perturbative 
expansion for which innovative optimization pro- 
cedures have been suggested. Among these are 
contour-improved (resummed) fixed-order per- 
turbation theory [14,15], effective charge and min- 
imal sensitivity schemes [16,17], the large- (3 ex- 
pansion [18,19], and combinations of these ap- 
proaches. They mainly distinguish themselves in 
how they deal with the fact that the perturbative 
series is truncated at an order where the missing 
part is not expected to be small. 

One could wonder how r decays may at all 
allow us to learn something about perturbative 
QCD. The hadronic decay of the r is dominated 
by resonant single particle final states. The corre- 
sponding QCD interactions that bind the quarks 
and gluons into these hadrons necessarily involve 
long distance scales, which are outside the domain 
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of perturbation theory. Indeed, it is the inclu- 
sive character of the sum of all hadronic r decays 
that allows us to probe fundamental short dis- 
tance physics. Inclusive observables like the total 
hadronic r decay rate R T can be accurately pre- 
dicted as function of a s (m 2 ) using perturbative 
QCD, and including small nonperturbative con- 
tributions within the framework of the Operator 
Product Expansion (OPE) [20]. In effect, R T is a 
doubly inclusive observable since it is the result 
of a summation over all hadronic final states at a 
given invariant mass and further over all masses 
between and m T . The scale m T lies in a com- 
promise region where a s (m 2 ) is large enough so 
that R T is sensitive to its value, yet still small 
enough so that the perturbative expansion con- 
verges safely and nonperturbative power correc- 
tions are small. 

If strong and electroweak radiative corrections 
are neglected, the theoretical parton level predic- 
tion for SUc{N c ), N c = 3 reads 

Rr = N c (\V ud \ 2 + \V US \ 2 ) = 3 , (7) 

and we can estimate a perturbative correction to 
this value of approximately 21%. One realizes the 
increase in sensitivity to a s compared to the Z 
hadronic width, where because of the three times 
smaller a s (M§ ) the perturbative QCD correction 
reaches only about 4%. 

The nonstrange inclusive observable R T can be 
theoretically separated into contributions from 
specific quark currents, namely vector (V) and 
axial-vector (A) ud and us quark currents. It is 
therefore appropriate to decompose 

Rr = Rr.V + Rr.A + Rr.S , (8) 

where for the strange hadronic width R T ,s vec- 
tor and axial-vector contributions are so far not 
separated because of the lack of the correspond- 
ing experimental information for the Cabibbo- 
supprcsscd modes. Parton-level and perturba- 
tive terms do not distinguish vector and axial- 
vector currents (for massless partons). Thus 
the corresponding predictions become R T via = 
(N c /2)\V ud \ 2 and R r , s = N C \V US \\ which add 
up to Eq. (7). 

A crucial issue of the QCD analysis at the r 
mass scale is the reliability of the theoretical de- 



scription, i.e., the use of the OPE to organize 
the perturbative and nonperturbative expansions, 
and the control of unknown higher-order terms 
in these series. A reasonable stability test is to 
continuously vary m T to lower values ^io < m T 
for both theoretical prediction and measurement, 
which is possible since the shapes of the r spec- 
tral functions are available. The kinematic factor 
that takes into account the r phase space sup- 
pression at masses near to m T is correspondingly 
modified so that ^/sq represents the new mass of 
the t. 

3.2. Theoretical prediction of R T 

According to Eq. (5) the absorptive parts of 
the vector and axial-vector two-point correla- 
tion functions n^ y , A (s), with the spin J of 
the hadronic system, are proportional to the 
r hadronic spectral functions with correspond- 
ing quantum numbers. The nonstrange ratio 
Rt,v+a can therefore be written as an integral of 
these spectral functions over the invariant mass- 
squared s of the final state hadrons [10] 

so 2 

Rt,v+a(sv) = 12tcS ew f—(l-±] x (9) 

o 

1 + 2— ImllW (s + ie) + Imn(°) (s + ie)\ , 

where iP 7 ) can be decomposed as n( J ) = 
\Vud\ 2 (n$v + II£2a) • The lower integration 
limit is zero because the pion pole is at zero mass 
in the chiral limit. 

The correlation function is analytic in the 
complex s plane everywhere except on the posi- 
tive real axis where singularities exist. Hence by 
Cauchy's theorem, the imaginary part of iP 17 ) is 
proportional to the discontinuity across the posi- 
tive real axis 

so 

J dsw(s)ImU(s) = ~ / dsw(s)U(s) , (10) 

|a|=s 

where w(s) is an arbitrary analytic function, 
and the contour integral runs counter-clockwise 
around the circle from s = s + ie to s = s — ie 
as indicated in Fig. 3. 
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Figure 3. Integration contour for the r.h.s. in 
Eq. (10). 

The energy scale sq = m 2 is large enough that 
contributions from nonpcrturbative effects are ex- 
pected to be subdominant and the use of the OPE 
is appropriate. The kinematic factor (1 — s/so) 2 
suppresses the contribution from the region near 
the positive real axis where Il( J ) (s) has a branch 
cut and the OPE validity is restricted due to large 
possible quark-hadron duality violations. 

The theoretical prediction of the vector and 
axial-vector ratio R t ,v/a can hence be written as 



R. 
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— -j\V u d\ 2 S-EW 
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with the massless perturbative contribution 5^°\ 
the residual non-logarithmic clcctroweak correc- 
tion 5' EW = 0.0010 [21], and the dimension D = 

2 perturbative contribution lrom quark 

masses. The term 8^ denotes the OPE con- 
tributions of mass dimension D 

. {O u d{p))v/A , s 
C ud ,v/A{s,n)— , (12) 



5 (d) _ y- 

J ud,V/A — / , 

dimO=D 



where the scale parameter \i separates the long- 
distance nonpcrturbative effects, absorbed into 



the vacuum expectation elements (O u d(^)), from 
the short-distance effects that are included in 
the Wilson coefficients C ud .v/A( s , L 1 )- Note that 
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3.2.1. The Perturbative Prediction 

The perturbative prediction used by the experi- 
ments follows the work of [14]. Effects from quark 
masses have been calculated in [22] and are found 
to be well below 1% for the light quarks. As a 
consequence, the contributions from vector and 
axial-vector currents coincide to any given order 
of perturbation theory and the results are flavor 
independent. 

For the evaluation of the perturbative series, 
it is convenient to introduce the analytic Adlcr 
function 



D(s) = -,s« 
as 



(13) 



The function D(s), calculated in perturbative 
QCD within the MS renormalization scheme, de- 
pends on a non-physical parameter fi occurring 
as ln(^ 2 /s). Furthermore it is a function of a s . 
On the other hand, since D(s) is connected to a 
physical quantity, the spectral function ImII(s), 
it cannot depend on the subjective choice of /z. 
This can be achieved if a s becomes a function of 
fi providing independence of D(s) of the choice 
of \x. Nevertheless, in the realistic case of a trun- 
cated series, some \x dependence remains and rep- 
resents an irreducible systematic uncertainty. 

The s dependence of the QCD coupling con- 
stant is obtained from the renormalization group 
equation (RGE) 



da s 
dins 



= (3(a s ) = -a 2 s J2^ 



(14) 



with a s = a s /n. Expressed in the MS renormal- 
ization scheme and for three active quark flavors 
at the t mass scale, the (3 n coefficients are known 
to four loops (n — 3) [23]. 

The perturbative expansion of the Adler func- 
tion can be inferred from the 3-loop calcula- 
tion of the e+e~ inclusive cross section ratio 
R e + e -(s) = cr(e + e~ — > hadrons (7))/er(e+e~ — > 
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/x+M" (7)) [24,25] 



1 OO 



(15) 



n=0 



where the K n (£) can be expressed [14] as function 
of K n and (3 n parameters (both known up to order 
n = 3), and £, an independent scale parameter. 
This leads to the perturbative expansion 



n=l 

with the functions 

*"'<«•> = h f T 



(16) 



(17) 



so 



3.2.2. Fixed-order perturbation theory 
(FOPT) 

Inserting the RGE solution for a s (s) into 
Eq. (16) to evaluate the contour integral, and col- 
lecting the terms with equal powers in a s leads to 
the familiar expression [14] 



5 (0) 



«s(£so) 



(18) 



where the g n are functions of K m<n and /3 m <n-i, 
and of elementary integrals with logarithms of 
power m < n in the integrand. Setting £ = 1 
and replacing all known and Ki coefficients by 
their numerical values, Eq. (18) simplifies to 

<5(°) = a s (s ) + (1.6398 + 3.5625) a 2 s (s ) 

+ (6.371 + 19.995) a 3 s (s ) (19) 

+ {Ka + 78.003) at(s ) 

+ (K 5 + 14.250 K 4 - 391.54) a 5 s (s ) , 

where for the purpose of systematic studies we 
have kept terms up to fifth order. When only two 
numbers are given in the parentheses, the first 
number corresponds to K n , and the second to g n . 

The FOPT series is truncated at given order de- 
spite the fact that parts of the higher coefficients 
g n >i{£,) are known to all orders and could be re- 
summed. These known parts are the higher (up 



to infinite) order logarithmic power terms of the 
expansion that are functions of the known /3„<3 
and K n < 3 only. In effect, beyond the use of the 
perturbative expansion of the Adler function (15), 
two approximations have been used to obtain the 
FOPT series (19): (i) the RGE (14) has been 
Taylor-expanded and terms higher than the given 
FOPT order have been truncated, and (ii) this 
Taylor expansion is used to predict a s (— s) on the 
entire \s\ — so contour. 

3.2.3. Contour-improved fixed-order per- 
turbation theory (CIPT) 

A more promising approach to the solution of 
the contour integrals Eq.(17) is to perform a di- 
rect numerical evaluation by means of single-step 
integration and using the solution of the RGE 
to four loops as input for the running a s (— £s) 
at each integration step [15,14]. It implicitly 
provides a partial resummation of the (known) 
higher order logarithmic integrals and improves 
the convergence of the perturbative series. While 
for instance the third order term in the expan- 
sion (19) contributes with 17% to the total (trun- 
cated) perturbative prediction, the correspond- 
ing term of the numerical solution amounts to 
only 6% (assuming a s (m 2 ) = 0.35). This nu- 
merical solution of Eq. (16) is referred to as 
contour-improved fixed-order perturbation theory 
(CIPT) in the following. Single-step integration 
also avoids the Taylor approximation of the RGE 
on the entire contour, since a s is iteratively com- 
puted from the previous step using the full known 
RGE. 

3.2.4. Other schemes and the value of K4 

Other approaches for evaluating the perturba- 
tive prediction have been presented, such as the 
effective charge perturbation theory (ECPT, see 
for instance Refs. [16,26]) and the large-/?o expan- 
sion [18,19]. Whereas these methods are of con- 
siderable theoretical interest, they are not suited 
for precision analyses [2]. 

ECPT has been largely used to estimate the 
value of the first uncalculated term K 4 [27]. 
Another approach to estimating K4 [28] is, within 
CIPT, to enhance the sensitivity to higher order 
perturbative terms by reducing the renormaliza- 
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tion scale £. Both methods yield K 4 ~ 27, but 
it has been shown [2] that the precision of this 
estimate is seriously limited (~ 100 %) by the 
lack of knowledge of the unknown higher order 
parameters in the perturbative series (K 5} . . .). 

Significant efforts are underway with the goal 
to calculate the K\ coefficient. Although the 
large number of five-loop diagrams that are 
needed to calculate the two-point current corre- 
lator at this order may appear discouraging, the 
results on two gauge invariant subsets are already 
available. The subset of order C(o^n^) was eval- 
uated long ago through the summation of renor- 
malon chains [29], while the much harder subset 
O(a^n'j) was recently calculated [30]. Following 
these investigations, the value K4 — 25 ± 25 is 
used in our analysis. 

3.2.5. Comparison of the perturbative 
methods 

To study the convergence of the perturbative 
series, we give in Tabic 1 the contributions of 
the different orders in PT to 5^ for the various 
approaches using a s (m 2 ) = 0.35. A geometric 
growth, K n <~ Kf l _ l / K n -2 1 is assumed for all 
unknown PT and RGE coefficients. In the case 
of CIPT the results are given for the various tech- 
niques used to evolve a s (s). 

Faster convergence is observed for CIPT com- 
pared to FOPT yielding a significantly smaller er- 
ror associated with the renormalization scale am- 
biguity. Our coarse extrapolation of the higher 
order coefficients could indicate that minimal sen- 
sitivity is reached at n ~ 5 for FOPT, while the 
series further converges for CIPT. Although the 
Taylor expansion in the CIPT integral exhibits 
significant deviations from the exact solution on 
the integration circle, the actual numerical effect 
from this on <j(°) is small (c/. second and third col- 
umn in Table 1). The convergence of the ECPT 
series is much worse than for FOPT and CIPT. 
Consequently, the difference between truncation 
at n — 4 and n = 6 may be significant. A similar 
instability may occur for the large-/3o expansion. 

The CIPT series is found to be better behaved 
than FOPT (as well as ECPT) and is therefore 
to be preferred for the numerical analysis of the 
r hadronic width. As a matter of fact, the differ- 



ence in the result observed when using a Taylor 
expansion and when truncating the perturbative 
series after integrating along the contour (FOPT) 
with the exact result at given order (CIPT) ex- 
hibits the incompleteness of the perturbative se- 
ries. However, it is even worse than that since 
large known coefficients are neglected in FOPT so 
that the difference between CIPT and FOPT may 
actually overstate the perturbative truncation un- 
certainty (certainly it is not a good measure of the 
latter uncertainty). This can be verified by study- 
ing the behavior of this difference for the various 
orders in perturbation theory given in Table 1. 
The CIPT-vs.-FOPT discrepancy increases with 
the addition of each order, up to order four where 
a maximum is reached. Adding the fifth order 
does not reduce the effect, and only beyond fifth 
order the two evaluations may become asymptotic 
to each other. As a consequence varying the un- 
known higher order coefficients and using the dif- 
ference between FOPT and CIPT as indicator of 
the theoretical uncertainties overemphasizes the 
truncation effect. 

3.3. Results 

It was shown in [11] that one can exploit 
the shape of the spectral functions to ob- 
tain additional constraints on a s (so) and — more 
importantly — on the nonperturbative effective 
operators. The r spectral moments at s = 
arc defined by 

m 2 ( 


where R®° v / A = Rt.v/a- The factor (1 - s/m^) k 
suppresses the integrand at the crossing of the 
positive real axis where the validity of the OPE 
less certain and the experimental accuracy is 
statistically limited. Its counterpart (s/m 2 ) 1 
projects upon higher energies. The spectral in- 
formation is used to fit simultaneously a s (m 2 ) 
and the effective operators (a s GG), pa s (qq) 2 and 
(On) for dimension D = 4, 6 and 8, respectively. 
Due to the large correlations between the differ- 
ently weighted spectral integrals, only five mo- 
ments are used as input to the fit. 
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Table 1 

Massless perturbative contribution to R r (m%) for the various methods considered, and at orders n > 1 
with a s (m 2 ) = 0.35. The value of K4 is set to 25, while all unknown higher order K n> 4 and (3 n >z 
coefficients are assumed to follow a geometric growth. Details are given in Ref. [2]. 

1° 



Pert. Method 


n = 1 


n = 2 


n = 3 


(n = 4) 


(n = 5) 


(n = 6) 




T 6 


FOPT (£ = f ) 


0.1114 


0.0646 


0.0365 


0.0159 


0.0010 


-0.0086 


0.2283 


0.2208 


CIPT (Taylor RGE, £ = 1) 


0.1573 


0.0317 


0.0126 


0.0042 


0.0011 


0.0001 


0.2058 


0.2070 


CIPT (full RGE, £ = 1) 


0.1524 


0.0311 


0.0129 


0.0046 


0.0013 


0.0002 


0.2009 


0.2025 


CIPT (full RGE, £ = 0.4) 


0.2166 


-0.0133 


0.0006 


-0.0007 


0.0010 


-0.0007 


0.2032 


0.2048 


ECPT 


0.1442 


0.2187 


-0.1195 


-0.0344 


-0.0160 


-0.0120 


0.2090 


0.1810 


Large-/?o expansion 


0.1114 


0.0635 


0.0398 


0.0241 


0.0155 


0.0093 


0.2388 


0.2636 



In analogy to R T , the contributions to the 
moments originating from perturbative and non- 
perturbative QCD are decomposed through the 
OPE. 

3.3.1. The ALEPH determination of 
c* s (m?) and nonperturbative con- 
tributions 

Combined fits to experimental spectral mo- 
ments and the extraction of a s (m^) together with 
the leading nonperturbative operators have been 
performed by ALEPH, CLEO and OPAL [4,5,7, 
12,13] using similar strategies and inputs. This 
analysis uses the final and complete data on 
branching fractions and spectral functions from 
ALEPH [5], yielding 

R T .v +A = 3.482 ± 0.014 , (21) 
R T ,v = 1-787 ± 0.011 ± 0.007 , (22) 
R t> a = 1-695 ± 0.011 ± 0.007 , (23) 

R t ,v-a = 0.092 ± 0.018 ± 0.014 , (24) 

where the second error originates from the V/A 
separation in final states with a KK pair, fully 
anticorrelated beetwen R r y and R t ,a- To re- 
duce the model dependence of the analysis, one 
fits simultaneously the nonperturbative opera- 
tors, which is possible since the correlations be- 
tween these and a s turn out to be small enough. 
The main theoretical uncertainties are due to K4 
(25±25) and to the renormalization scale, which 
is varied around m T from 1.1 to 2.5 GeV (the vari- 
ation over half of the range taken as systematic 



uncertainty). 

The fit results are given in Table 2. There is a 
remarkable agreement within statistical errors be- 
tween the a s (m 2 ) determinations using the vec- 
tor and axial-vector data. This provides an im- 
portant consistency check of the results, since the 
two corresponding spectral functions are experi- 
mentally independent and manifest a quite dif- 
ferent resonant behavior. However it must be 
mentioned that the a s (m 2 ) determination using 
either the V and A spectral functions is more 
dependent on the validity of the OPE approach 
since their nonperturbative contributions are sig- 
nificantlty larger than for V + A. Indeed the lead- 
ing nonperturbative contributions of dimension 
D = 6 and D = 8 approximately cancel in the 
inclusive sum. This cancellation of the nonper- 
turbative terms increases the confidence in the 
a s (m 2 ) determination from the inclusive V + A 
observables. Averaging CIPT and FOPT, the re- 
sult quoted by ALEPH is 

a s (m 2 T ) = 0.340 ± 0.005 cxp ± 0.014 th , (25) 

The gluon condensate is determined by the first 
k = 1, I = 0, 1 moments, which receive lowest 
order contributions. The values obtained in the 
V and A fits are not very consistent, which could 
indicate problems in the validity of the OPE ap- 
proach used once the nonperturbative terms be- 
come significant. Taking the value obtained in the 
V+A fit, where nonperturbative effects are small, 
and adding as systematic uncertainties half of the 
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Table 2 

Results [5] for a s (m^) and the nonperturbative contributions for vector, axial-vector and V+A combined 
fits using the corresponding experimental spectral moments as input parameters. Where two errors are 
given the first is experimental and the second theoretical. The 5^ term is theoretical only with quark 
masses varying within their allowed ranges (see Ref. [2]). The quark condensates in the 5^ term are 
obtained from PCAC, while the gluon condensate is determined by the fit. The total nonperturbative 
contribution is the sum (5np = + + . Full results are listed only for the CIPT perturbative 
prescription, except for a s (m^) where the results using both CIPT and FOPT are given (See Ref. [1] 
for the complete results). 



Parameter 


Vector (V) 


Axial- Vector (A) 


V + A 


a s (m 2 T ) (CIPT) 
a s (m?) (FOPT) 


0.355 ±0.008 ±0.009 
0.331 ±0.006 ±0.012 


0.333 ± 0.009 ± 0.009 
0.327 ±0.007 ±0.012 


0.350 ± 0.005 ± 0.009 
0.331 ±0.004 ±0.012 


5 (2) (CIPT) 


(-3.3 ±3.0) x 10" 4 


(-5.1 ±3.0) x 10" 4 


(-4.4 ±2.0) x 10~ 4 


{a 3 GG) (GeV 4 ) (CIPT) 


(0.4 ±0.3) x 10 -2 


(-1.3 ±0.4) x 10~ 2 


(-0.5 ±0.3) x 10~ 2 


S W (CIPT) 


(4.1 ± 1.2) x 10~ 4 


(-5.7 ±0.1) x 10~ 3 


(-2.7 ±0.1) x 10 -3 


5 (6) (CIPT) 


(2.85 ±0.22) x 10~ 2 


(-3.23 ±0.26) x 10~ 2 


(-2.1 ± 2.2) x 10~ 3 


5 (8) (CIPT) 


(-9.0 ±0.5) x 10~ 3 


(8.9 ±0.6) x 10~ 3 


(-0.3 ±4.8) x 10~ 4 


Total 5 NP (CIPT) 


(1.99 ±0.27) x 10~ 2 


(-2.91 ±0.20) x 10~ 2 


(-4.8 ± 1.7) x 10~ 3 



X 2 /DF (CIPT) 0.52 4.97 3.66 



difference between the vector and axial- vector fits 
as well as between the CIPT and FOPT results, 
ALEPH measures the gluon condensate to be 

(a s GG) = (0.001 ± 0.012) GeV 4 . (26) 

This result docs not provide evidence for a 
nonzero gluon condensate, but it is consistent 
with and has comparable accuracy to the inde- 
pendent value obtained using charmonium sum 
rules and e + e~ data in the charm region, (0.011± 
0.009) GeV 4 in a combined determination with 
the c quark mass [31]. 

The approximate cancellation of the nonper- 
turbative contributions in the V + A case was 
predicted [10] for D — 6 assuming vacuum satu- 
ration for the matrix elements of four-quark oper- 
ators, which yields 5 ( y ] /5f = -7/11 = -0.64, in 
fair agreement with the result — 0.90± 0.18. The 
estimate [10] for c^ 6) = (2.5 ± 1.3) x 10~ 2 agrees 
with the experimental result. 

The total nonperturbative V + A correction, 
<5np,v+a = (-4.3 ± 1.9) x 10~ 3 , is an order of 



magnitude smaller than the corresponding values 
in the V and A components, Jnp.v = (2.0 ±0.3) x 
10~ 2 and <5 NPjA = (-2.8 ±0.3) x 10~ 2 . 

3.3.2. Running of a s {s) below 

Using the spectral functions, one can simulate 
the physics of a hypothetical r lepton with a mass 
y/so smaller than m T through Eq. (9). Assum- 
ing quark-hadron duality, the evolution of R T (so) 
provides a direct test of the running of a s (so), 
governed by the RGE /3-function. On the other 
hand, it is also a test of the stability of the OPE 
approach at small scales. The studies performed 
in this section employ only CIPT. 

The measured function R T y + A(so) is plotted 
in the left hand plot of Fig. 4 together with the 
theoretical prediction using the results of Table 2. 
The correlations between two adjacent points in 
so are large as the only new information is pro- 
vided by the small mass difference between the 
two points and the slightly modified weight func- 
tions under the integral. Moreover the correla- 
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0.5 0.75 1 1.25 1.5 1.75 2 2.25 2.5 2.75 3 1 1.25 1.5 1.75 2 2.25 2.5 2.75 3 



s (GeV 2 ) s (GeV 2 ) 

Figure 4. Left: The ratio R T y+A versus the square "r mass" so- The curves are plotted as error bands 
to emphasize their strong point-to-point correlations in sq. Also shown is the theoretical prediction using 
CIPT and the results for R T y + A and the nonperturbative terms from Table 2. Right: The running of 
a s (so) obtained from the fit of the theoretical prediction to i? Tj y+A(so) using CIPT. The shaded band 
shows the data including only experimental errors. The curve gives the expected four-loop RGE evolution 
for three flavors. 



tions are reinforced by the original experimental 
and theoretical correlations. Below 1 GeV 2 the 
error of the theoretical prediction of R T y+A(so) 
starts to blow up because of the increasing sensi- 
tivity to the unknown K4 perturbative term; er- 
rors of the nonperturbative contributions are not 
contained in the theoretical error band. Figure 4 
(right) shows the plot corresponding to Fig. 4 
(left), translated into the running of a s (so). Only 
experimental errors are shown. Also plotted is 
the four-loop RGE evolution using three quark 
flavors. 

It is remarkable that the theoretical predic- 
tion using the parameters determined at the r 
mass and R t ,v+a{sq) extracted from the mea- 
sured V + A spectral function agree down to 
s ~ 0.8 GeV 2 . The agreement is good to about 
2% at 1 GeV 2 . This result, even more directly il- 
lustrated by the right hand plot of Fig. 4, demon- 
strates the validity of the perturbative approach 
down to masses around 1 GeV, well below the r 
mass scale. The agreement with the expected 
scale evolution between 1 and 1.8 GeV is an inter- 



esting result, considering the relatively low mass 
range, where a s is seen to decrease by a factor 
of 1.6 and reaches rather large values ~ 0.55 at 
the lowest masses. This behavior provides confi- 
dence that the a s (m 2 ) measurement is on solid 
phenomenological ground. 

3.3.3. Final assessment on the a s (m^) de- 
termination 

Although this evaluation of a s (m 2 ) represents 
the state-of-the art, several remarks can be made: 

• The analysis is based on the ALEPH spec- 
tral functions and branching fractions, en- 
suring a good consistency between all the 
observables, but not exploiting the full ex- 
perimental information currently available 
from other experiments. Because the re- 
sult on a s (m 2 ) is limited by theoretical un- 
certainties, one should expect only a small 
improvement of the final error in this way, 
however it can influence the central value. 

• One example for this is the evaluation of 
the strange component. Some discrepancy 
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is observed between the ALEPH measure- 
ment of the (Kitit)~v mode and the CLEO 
and OPAL results. Although this could 
still be the result of a statistical fluctuation, 
their average provides a significant shift in 
the central value compared to using the 
ALEPH number alone. Another improve- 
ment is the substitution of the measured 
branching fraction for the K~v mode by the 
more precise value predicted from r-fi uni- 
versality. Both operations have the effect 
to increase R T ,s, the ratio of the r decay 
width into strange hadronic final states to 
the electronic width, from 0.1603 ± 0.0064, 
as obtained by ALEPH, to 0.1686 ± 0.0047 
for the world average. 

• One can likewise substitute the world av- 
erage value for the universality-improved 
value of the electronic branching fraction, 
B™ = (17.818 ± 0.032)%, to the ALEPH 
result, B e = (17.810 ± 0.039)%, with little 
change in the central value, but some im- 
provement in the precision. 

From this analysis, one finds the new value for 
the nonstrange ratio, 

Rt,v+a = Rr — Rt.S 

= (3.640 ± 0.010) - (0.1686 ±0.0047) 
= 3.471 ±0.011. (27) 

The result (27) translates into the following de- 
termination of a s (m 2 ) from the inclusive V + A 
component using the CIPT approach 

a s {m 2 T ) = 0.345 ± 0.004 oxp ± 0.009 th , (28) 

with improved experimental and theoretical pre- 
cision over the ALEPH result. Most of the the- 
oretical uncertainty originates from the limited 
knowledge of the perturbative expansion, only 
predicted to third order. Following the dicussion 
above we take the result from the CIPT expan- 
sion, not introducing any additional uncertainty 
spanning the difference between FOPT and CIPT 
results. The dominant theoretical errors are from 
the uncertainty on K4 and from the renormaliza- 
tion scale dependence, both covering the effect of 
truncating the series after the estimated fourth 
order. 



3.3.4. Evolution to Mf 

It is customary to compare a s values, obtained 
at different renormalization scales, at the scale of 
the Z-boson mass. 

The evolution of the a s (m%) measurement 
from the inclusive V + A observables given in 
Eq. (28), based on Runge-Kutta integration of 
the RGE (14) to N 3 LO , and three-loop quark- 
flavor matching, gives 

a 3 [M%) = 0.1215 (4 cxp ) (10) th (5) evol , 

= 0.1215 ± 0.0012 . (29) 

The first two errors originate from the a s (m 2 ) de- 
termination given in Eq. (28). The last error re- 
ceives contributions from the uncertainties in the 
c-quark mass (0.00020, m c varied by ±0.1 GeV) 
and the 6-quark mass (0.00005, mj varied by 
±0.1 GeV), the matching scale (0.00023, (j, var- 
ied between 0.7 m q and 3.0 m q ), the three-loop 
truncation in the matching expansion (0.00026) 
and the four-loop truncation in the RGE equation 
(0.00031), where we used for the last two errors 
the size of the highest known perturbative term 
as systematic uncertainty. These errors have been 
added in quadrature. The result (29) is a deter- 
mination of the strong coupling at the Z mass 
scale with a precision of 1%. 

The evolution path of a s {m%.) is shown in the 
upper plot of Fig. 5. The two discontinuities are 
due to the quark-flavor matching at fj, = 2m q . 
One could prefer to have an (almost) smooth 
matching by choosing /i = m q . However, in this 
case, one must first evolve from m T down to m c 
to match the c-quark flavor, before evolving to 
nib. The effect on a s (M§) from this ambiguity 
is within the assigned systematic uncertainty for 
the evolution. 

The comparison with the other determinations 
of a s (M§) is given in Fig. 5 using compiled re- 
sults from Ref. [32]. 

3.3.5. A measure of asymptotic freedom 
between and M% 

The r-decay and Z-width determinations have 
comparable accuracies, which are however very 
different in nature. The r value is dominated by 
theoretical uncertainties, whereas the determina- 
tion at the Z resonance, benefiting from the much 
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DIS (e/|Xf F 2 ) 

[1 .9-1 5.2] 



a 



10 

(i scale (GeV) 



Figure 5. Top : The evolution of a s {m 2 T ) (28) to 
higher scales \i using the four-loop RGE and the 
3-loop matching conditions applied at the heavy 
quark-pair thresholds (hence the discontinuities 
at 2m c and 2m(,). The evolution is compared with 
other independent measurements (see text) cover- 
ing scales varying over more than two orders mag- 
nitude. The experimental values are taken from 
the compilation [32]. Bottom : The corresponding 
extrapolated a s values at Mz- The shaded band 
displays the r decay result within errors. 

larger energy scale and the correspondingly small 
uncertainties from the truncated perturbative ex- 
pansion, is limited by the experimental precision 
on the electroweak observables, essentially the ra- 
tio of leptonic to hadronic peak cross sections. 
The consistency between the two results provides 
the most powerful present test of the evolution of 
the strong interaction coupling, as it is predicted 
by the nonabelian nature of the QCD gauge the- 
ory. This test extends over a range of s spanning 
more than three orders of magnitude. The differ- 
ence between the extrapolated r-decay value and 
the measurement at the Z is: 



which agrees with zero with a relative precision 
of 2.4%. 

In fact, the comparison of these two values is 
valuable since they are among the most precise 
single measurements and they are widely spaced 
in energy scale. Thus it allows one to perform 
an accurate test of asymptotic freedom. Let us 
consider the following evolution estimator [33] for 
the inverse of a s (s), 

a- 1 (si) - aj 1 (s 2 ) 



r(si,s 2 ) 



(31) 



lnsi — lns2 

which reduces to the logarithmic derivative of 
a~ 1 (s) when si — ► S2, 



2tt0{s) 



dim 



2/3o 

7T 



01 Ota 

1 + 7T — 

P0 7T 



(32) 



a T s {M%) - af (Aff ) = 0.0029 (10) T {27) z 



(30) 



with the notations of Eq. (14). At first order, the 
logarithmic derivative is driven by 0q. 

The r and Z experimental determinations of 
a s {s) yield the value 

r exp {m 2 T , Mf ) = 1.405 ± 0.053 , (33) 

which agrees with the prediction using the RGE 
to N 3 LO, and 3-loop quark- flavor matching, 

r QC B{m 2 T , M 2 Z ) = 1.353 ± 0.006 . (34) 

To our knowledge this is the most precise exper- 
imental test of the asymptotic freedom property 
of QCD at present. It can be compared to an in- 
dependent determination [33], using an analysis 
of event shape observables at LEP between the 
Z energy and 207 GeV, r(Mf , (207 GeV) 2 ) = 
1.11 ± 0.21, for a QCD expectation of 1.27. 

4. Conclusions 

Using the ALEPH spectral functions and 
branching ratios, complemented by other avail- 
able measurements, and a revisited analysis of 
the theoretical framework, the value a s {m 2 ) = 
0.345 ± 0.004 cxp ± 0.009th is obtained. Taken to- 
gether with the determination of a s (M§) from 
the global electroweak fit, this result leads to the 
most accurate test of asymptotic freedom: the 
value of the logarithmic slope of aj 1 (s) is found 
to agree with QCD at a precision of 4%. 
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